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Abstract. We review the mechanisms of low-temperature electron transport across a quantum dot 
weakly coupled to two conducting leads. Conduction in this case is controlled by the interaction 
between electrons. At temperatures moderately lower than the single-electron charging energy of 
the dot, the linear conductance is suppressed by the Coulomb blockade. Upon further lowering of 
the temperature, however, the conductance may start to increase again due to the Kondo effect. 
This increase occurs only if the dot has a non-zero spin 5*. We concentrate on the simplest case 
of S = 1/2, and discuss the conductance across the dot in a broad temperature range, which 
includes the Kondo temperature. Temperature dependence of the linear conductance in the Kondo 
regime is discussed in detail. We also consider a simple (but realistic) limit in which the differential 
conductance at a finite bias can be fully investigated. 



1. Introduction 

In this review, we define a quantum dot as a small paddle of electron liquid which 
is attached by tunnel junctions to two massive conductors (leads). The number of 
electrons on the dot in equilibrium N can be controlled by varying the voltage 
applied to a gate - an auxiliary electrode coupled to the dot only capacitively 
(no tunneling occurs between the dot and gate). If the conductances of the lead- 
dot junctions are small compared to the conductance quantum /h, and if the 
temperature T is low enough, then N is an integer at almost any gate voltage Vg. 
Exceptions are small intervals around the discrete set of values of V^, at which 
an addition of a single electron to the dot does not change the electrostatic en- 
ergy. Such a degeneracy between different charge states of a quantum dot allows 
for an activationless electron transfer through it, whereas for all other values of 
Vg the activation energy for the conductance G across the dot is finite. The re- 
sulting oscillatory dependence G{Vg) is the hallmark of the Coulomb blockade 
phenomenon |1|. The contrast between the low- and high-conductance regions 
(Coulomb blockade valleys and peaks, respectively) gets sharper at lower tem- 
peratures. The described pattern of the G{Vg,T) dependence is clearly observed 
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down to the lowest attainable temperatures in experiments on tunneling through 
small metallic islands |2| However, small quantum dots formed in GaAs het- 
erostructures display a different behavior |3]| : in some Coulomb blockade valleys 
the dependence G{T) is not monotonic and has a minimum at a finite temperature. 
This minimum is similar in origin |4 1 to the well-known non-monotonic tempera- 
ture dependence of the resistivity of a metal containing magnetic impurities |5l - 
the Kondo effect. 

We present here an introduction to the theory of Kondo effect in quantum dots, 
concentrating on the so-called Constant Interaction Model. Despite its simplicity, 
the model allows one to describe the main conduction mechanisms, including the 
most important version of the Kondo effect in quantum dots. 

2. Constant Interaction Model 

The Hamiltonian of interacting electrons confined to a quantum dot has the fol- 
lowing general form: 



Here an operator d/-^ creates an electron with spin s in the orbital state </>j(r); 
Hij = is an Hermitian matrix describing the single-particle spectrum, and 
the matrix elements Hijki depend on the electron-electron interaction potential 
V{ri — and on the chosen basis of orbital states: 



Further simplification of the Hamiltonian Q is possible if the following condi- 
tions are met: 

(/) the electron-electron interaction within the dot is not too strong (the gas pa- 
rameter characterizing the electron-electron interaction in a non-ideal Fermi gas 
rs < 1) so that the Fermi liquid description is applicable; 

(//) the quasiparticle spectrum is not degenerate near the Fermi energy (which is 
satisfied, in general, if the motion of electrons within the dot is chaotic); 
(///) the dot is in the metallic regime, i.e. the Thouless energy of the dot Et 
exceeds the mean single-particle level spacing 6E. 

The Thouless energy Et is well-defined if the electron motion within it is 
chaotic due to either disorder or irregular shape of the dot. For a dot of linear size 
L it equals 




(1) 





1->L 



(3) 
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where I if the elastic mean free path and D is the diffusion constant in a disordered 
dot. The average level spacing 6E can be estimated as 

6E ~ l/i^dL"^. (4) 

Here d is the dimensionality of the system (d = 2 for a quantum dot formed in a 
two-dimensional electron gas in a semiconductor heterostructure, and d = 3 for 
a metallic nanoparticle); Ud is the density of states of the material. For a ballistic 
electron motion (/ ^ L) the ratio 

Et/5E ~ L/Aj7 ~ iVi/'^ 

is large for large number of electrons N on the dot. 

If the conditions are met, then the Random Matrix Theory (RMT) is 

a good starting point for describing non-interacting quasiparticles, see ||6l for a 
review; the matrix elements Tiij belong to a Gaussian ensemble fT). The matrix 
elements do not depend on spin, and therefore each eigenvalue of the matrix 
Tiij represents a spin-degenerate energy level. The spacings le^+i — e^l between 
consecutive levels obey the Wigner-Dyson distribution |7l; the average value of 
|e„+i - e„| is 5E. 

It turns out jH that the vast majority of the matrix elements Tiijki are small. 
In the Umit Et/SE oo, only relatively few "most diagonal" elements re- 
main finite. In this limit the interaction part Hint of the Hamiltonian of a non- 
superconducting dot may be cast in the form \9]: 

H,nt = Ec(N -N)^ -EsS\ (5) 

Here 

N = Y,disdns, ^ = T.^ns^<s' (6) 

ns nss' 

are the operators of total number of electrons and total spin on the dot. Obviously, 
the Hamiltonian ^ is invariant with respect to a change of the single-particle 
basis (j)i{r). 

The first term in Eq. ^ represents the electrostatic energy. In the conventional 
equivalent circuit picture, see Fig. [0 the charging energy Ec is related to the 
total capacitance C of the dot, Ec = e^/2C, and the dimensionless parameter 
M is proportional to the gate voltage, J\f = CgVg/e, where Cg is the capacitance 
between the dot and the gate. For a mesoscopic (Ai? <C L) conductor, the charging 
energy is large compared to the level spacing 5E. Indeed, using the estimate C ~ 
kL, where k is the dielectric constant, we find 

5E nhvp (^i?) s ! ( ) 
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where vpis the Fermi velocity, and is the conventional gas parameter. Except an 
exotic case of an extremely weak interaction, having a large number of electrons 
N ~ {L/\fY ^ lona quantum dot guarantees that the inequality Ec 2> 6E 
is satisfied in dimensions d = 2,3. For the smallest quantum dots formed in 
GaAs heterostructures, Ec/6E ~ 10 |3|; for a metallic nanoparticle with L ~ 
5 nm this ratio is about 100, see, e.g., |10|. We should mention also that in the 
d = 1 case, the estimate Q is of little help, as even a weak electron-electron 
interaction results in a formation of a Luttinger liquid 1 11 1, significantly affecting 
the nomenclature of the elementary excitations. In the case of a single-mode finite- 
length Luttinger liquid, the elementary excitations can be viewed as Id plasmon 
waves in a confined geometry, so the corresponding level spacing is 5E ~ Ec- 
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Figure 1. Equivalent circuit of a dot connected to two leads by tunnel junctions (with conductances 
Gl and Gr) and coupled via a capacitor (with capacitance Cg) to the gate. The total capacitance 
of the dot C = Cl + Cr + Cg. 

The second term in Eq. (|5Jl is the exchange interaction characterized by the 
exchange constant Es- This constant is small in the case of weak electron-electron 
interaction: Es ~ rsSE (there is an additional ln(l/rs) factor in this estimate in 
the case of a 2d dot). Smallness of the ratio Es/6E guarantees the absence of a 
macroscopic (proportional to the volume of the dot) spin in the ground state, in 
accordance with the well-known Stoner criterion for the itinerant magnetism lIT^ . 
If all the orbital levels were equidistant, then the spin of an even-electron state 
would be zero, while an odd-electron state would have spin 1/2. However, the 
level spacings |e„+i — e^l are random. If the spacing between some orbital levels 
is accidentally small, the dot may acquire a spin fOll exceeding 1/2. At small r^, 
however, such states are rare, and in the following we disregard this possibility. 
This allows us to neglect the exchange term in the interaction Hamiltonian 

Keeping in mind that Eq. (|5Jl is invariant with respect to the rotations of the 
basis of single-particle states, we can pick up the basis in which the first term in 
the Hamiltonian Q is diagonal. Eq. (ISj then reduces to the Hamiltonian of the 
Constant Interaction Model, 

Hdot = Yl ^ndisdns + Ec{n-M)'' , (8) 

ns 
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which is widely used in the analysis of the experimental data. We stress that the 
single-particle level spacings |e„+i — e„| ~ dE here are small compared to the 
charging energy, and each single-particle state is spin-degenerate. 

Electron transport through the dot occurs via two dot-lead junctions, see Fig. [2 
As long as the conductances Gl and Gr of the left (L) and right (R) junctions 
are small compared to e^//i, one may use the tunneling Hamiltonian to describe 
the junctions: 

E WLA. + H.c., a = R,L. (9) 

akns 

Here cj^,^^ are the creation operators of electrons in the leads, and takn are the 
tunneling matrix elements (tunneling amplitudes) connecting the state n in the 
dot with the state k in the lead a. Usually it is adequate to consider the leads as 
reservoirs of free electrons: 

'^a = E ^ak'^ks'^aks ^^^^ 
ks 

with continuous spectra ^^k- In the following we will characterize these spectra 
by the density of states v, same for both leads. 

The randomness of states n translates into the randomness of the tunneling 
amplitudes tatn- Similar to the single-particle spectrum of the dot, the tunneling 
amplitudes can be studied in the RMT framework. Regardless the details of the 
statistical properties of takn, the corresponding tunneling probabilities can be 
related to the average conductances Ga of the junctions: 

EJ=3^*^. (11) 

The tunneling amplitudes determine also the rate Fa for an electron occupying a 
given discrete level in the dot to escape into the lead a = L,R. The escape rates 
are determined by the tunneling amplitudes. With the help of Eq. (fTlTi the average 
values of the rates are estimated to be Tq, ~ {hGa/ e'^)5E. 

3. Rate equations and conductance across the dot 

At high temperatures, T ^ Ec, charging energy is negligible compared to the 
thermal energy of electrons. Therefore the conductance Goo in this regime is not 
affected by charging and, independently of the gate voltage, 

' ' (12) 



Goo Gl Gji 

At lower temperatures. 



5E <^ Ec, (13) 
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electron transport across the dot starts to depend on N. The conductance is not 
suppressed only within narrow regions - Coulomb blockade peaks, i.e., at the gate 
voltages tuned sufficiently close to the points of charge degeneracy, 

W-N*\<T/Ec-, (14) 

here TV* is a half-integer number. 

We will demonstrate this now using the method of rate equations. In addition 
to the constraint Eq. (fT3l . we will assume that the inelastic electron relaxation 
rate l/r^ within the dot is large compared to the escape rates F^. In other words, 
electron transitions between discrete levels in the dot occur before the electron es- 
capes to the leads^. Under this assumption the tunnelings across the two junctions 
can be treated independently of each other. The condition (fUTi . on the other hand, 
allows us to take into account only two charge states of the dot which are almost 
degenerate in the vicinity of a Coulomb blockade peak. For M close to J\f* these 
are the states with N = M* ±1/2 electrons on the dot. Hereafter we denote these 
states as |A^) and + 1). Finally, the condition (fT3l enables us to replace the 
discrete single-particle levels within the dot by a continuum with the density of 
states 1/6E. 

Applying the Fermi Golden Rule and using the described simplifications we 
may write the current from the lead a into the dot as 

la = -^'^\takn\'^^iCk + eVa + En - en - En+i) 
kns 

X {PNf{Ck)[l - /(en)] - Pjv+i/(en)[l - /fe)]} . 

Here f{uj) = [exp(a;/T) + 1]^^ is the Fermi function, Va is the potential on 
the lead a, see Fig.^ En and En+i are the electrostatic energies of the charge 
states |A^) and [A^ + 1), and P/v and P/v+i are the probabilities to find the dot 
in these states. Replacing the summations over n and k by integrations over the 
corresponding continua, we find 

Ia = — [PnF{Ei -Eo- i^a) - Pn+iF{Eq - Ei + li^)] (15) 
e 

with F{uj) = u;/[exp(Lij/T) — 1]. In equilibrium, the current Eq. (fTSl is zero 
by the detailed balance. When a finite current flows through the junction the 
probabilities P/v and Pn+i deviate from their equilibrium values. In the steady 
state, the currents across the two junctions satisfy 

I = Il = -Ir- (16) 

^ Note that a finite inelastic relaxation rate requires inclusion of mechanisms beyond the 
Constant Interaction Model, e.g., electron-phonon collisions. 
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From Eqs. (fT5t and (fT6l . supplemented by the obvious normalization condition 
Ptv + Pn+1 = 1. one finds the current across the dot / in response to the applied 
bias V = Vl — Vr. The linear conductance across the dot G = limy^o dl/dV is 
then given by 

G-G Ec{M-M*)/T 

Here a half-integer J\f = J\f* corresponds to the Coulomb blockade peak. In the 
Coulomb blockade valleys {M ^ N*), conductance falls off exponentially with 
the decreasing temperature, and all the valleys behave exactly the same way. 

According to Eq. (flTl . all the Coulomb blockade peaks have the same height 
G{M*) = Goo/ 2 and are symmetric about J\f = M*. Such equivalence of all 
the peaks occurs because for T 3> 6E, see Eq. (fT3b . many energy levels in the 
dot make a significant contribution to the conductance. Decrease of temperature 
below 6E brings about a modification of Eq. (I17t in the vicinity of the Coulomb 
blockade peaks \J\f — M* \ <C 6E/Ec- In this regime the main contribution to the 
conductance comes from tunneling through a single energy level in the dot. The 
dependence G{Af, T) can again be found from the proper rate equations with the 
result m 

0.-gJ4'J^. X^'-E^iM-N'). (18, 

In writing Eq. (fTSl we neglected mesoscopic fluctuations of the peaks heights fl4ll . 
Note that according to Eq. (fTSl the shape of the peak is asymmetric. The asymme- 
try comes about from the difference in the spin states of the dots in two adjacent 
Coulomb blockade valleys: Eq. (fTSl is written for a peak separating the valley 
with A'^ = TV* — 1/2 = odd number of electrons at X < from the valley 
with N = M* + 1/2 = even at X > 0. The validity of Eq. (CH) is restricted 
to the vicinity of the Coulomb blockade peak, \J\f — N*\ <C 5E/Ec, and for 
temperatures in the interval ^ T ^ 5E. When the temperature falls below 
the escape rates Fq the Coulomb blockade peaks are no longer well-defined due 
to the onset of the Kondo effect. 



4. Activationless transport through a blockaded quantum dot 

According to the rate equations theory |15|, at low temperatures, T <^ Eq, 
conduction through the dot in the Coulomb blockade valleys is exponentially sup- 
pressed. This suppression occurs because the process of electron transport through 
the dot involves a real transition to the state in which the charge of the dot differs 
by e from the thermodynamically most probable value. The probability of such 
fluctuation is proportional to exp[— £'c|AA — M*\/T\, which explains the con- 
ductance suppression, see Eq. (flTl . Going beyond the lowest-order perturbation 
theory in conductances Ga allows one to consider processes in which states of the 
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dot with a "wrong" charge participate in the tunneling process as virtual states. 
The existence of these higher-order contributions to the tunneling conductance 
was first envisioned by Giaever and Zeller fT6l . The first quantitative theory of 
this effect, however, was developed much later fV7]. 

The leading contributions to the activationless transport, according to Refs. 1 17 1, 
are provided by the processes of inelastic and elastic co-tunneling. Unlike the 
sequential tunneling, in the co-tunneling mechanism, the events of electron tun- 
neling from one of the leads into the dot, and tunneling from the dot to the other 
lead occur as a single quantum process. 



4.1. INELASTIC CO-TUNNELING 



In the inelastic co-tunneling mechanism, an electron tunnels from the lead into 
one of the vacant single-particle levels in the dot, while it is an electron occupying 
some other level that tunnels out of the dot, see Fig. I3a). As a result, transfer 
of charge e between the leads is accompanied by a simultaneous creation of an 
electron-hole pair in the dot. 




Figure 2. Examples of the co-tunneling processes. 

(a) inelastic co-tunneling: transferring of an electron between the leads leaves behind an elec- 
tron-hole pair in the dot; (b) elastic co-tunneling; (c) elastic co-tunneling with a flip of spin. 



Here we will estimate the contribution of the inelastic co-tunneling to the con- 
ductance deep in the Coulomb blockade valley, i.e., at almost integer J\f. Consider 
an electron that tunnels into the dot from the lead a. If the energy of the electron 
relative to the Fermi level lo is small compared to the charging energy, lu <g; Ec, 
then the energy of the virtual state involved in the co-tunneling process is close 
to Ec. The amplitude ylj„ of the inelastic transition via this virtual state is then 
given by 

f* f 

Ain akn a'k'n' ^^Q\ 

The initial state of this transition has an extra electron in the single-particle state 
k in the lead a, while the final state has an extra electron in the state k' in the lead 
a' and an electron-hole pair in the dot (state n is occupied, state n' is empty). 
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Given the energy of the initial state u, the number of available final states 
can be estimated from the phase space argument, familiar from the calculation 
of the quasiparticle lifetime in the Fermi liquid theory fTSll . For uj S> 5E this 
number is of the order of {uj/SE)"^. Since the typical value of uj is T, the inelastic 
co-tunneling contribution to the conductance can be estimated as 

Using Eqs. ( fT9l and (fTTT l we find I17il 



G^n- ^GlGr[ — ] . (20) 



h 

A comparison of Eq. d20b with the result of the rate equations theory Eq. (fTTb 
shows that the inelastic co-tunneling takes over the thermally-activated hopping 
at moderately low temperatures 

-1 



T<Tin = Ec 



In 



e'/h 



Gl + Gr 



(21) 



The smallest energy of the electron-hole pair is of the order of 5E. At temperatures 
below this threshold the inelastic co-tunneling contribution to the conductance 
becomes exponentially small. It turns out, however, that even at much higher 
temperatures this mechanism becomes less effective than the elastic co-tunneling. 



4.2. ELASTIC CO-TUNNELING 

In the process of elastic co-tunneUng, transfer of charge between the leads is not 
accompanied by the creation of an electron-hole pair in the dot. In other words, 
occupation numbers of single-particle energy levels in the dot in the initial and 
final states of the co-tunneling process are exactly the same, see Fig. I3b). Here 
we will estimate the elastic co-tunneling contribution to the conductance near the 
edge of a Coulomb-blockade valley, 

5E/ Ec -M* <.l/2. (22) 

Under this condition, the average number of electrons on the dot N = J\f* + 
1/2. A cost in electrostatic energy for a removal of a single electron from the 
dot, S_ = 2Ec{M — A/"*), see Eq. (jSJl, is small compared to the cost to add an 
electron Ej^ ?a 2Ec [here we took into account the second inequality in Eq. (l22l l. 
Therefore, only hole-like virtual states with — 1 electrons on the dot contribute 
to the co-tunneling amplitude: 

Ad _ Sr^ An An _ ^akn^a'k'n 
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Here represent amplitudes of the processes in which a hole is virtually created 
on nth single-particle level. Creation of a hole is only possible if the level is 
occupied, hence the restriction of the sum in Eq. (123 1 to the levels below the 
Fermi level (en < 0). 

The elastic co-tunneling contribution to the conductance is 



Gei = — 
TTh 



(24) 



If the dot-leads junctions are point contacts, then Gei exhibits strong mesoscopic 
fluctuations. Indeed, tunneling matrix elements entering Eq. depend on the 
values of the electron wave functions at the points Tq, of the contacts, takn oc 
¥'n(i'Q). In the RMT model, briefly discussed in Section|2l the electron wave func- 
tions in the dot ^pn{'^) are random and uncorrelated, ^n{^o)'^n'{^a') Snn'^aa' 



Therefore, the partial amplitudes vl^^ „ are random and uncorrelated as well: 



L^RJ 



A 



L^R 



Eqs. (I23t . d24l i and (fTTT l then yield for the average conductance 



G, 



el 



en<0 



L^R\ ~ g2 



—glGr 

en<0 



5E 



n 2 



E- + |e„| 



Since under the conditions d22l the number of terms making significant contri- 
bution to the sum over n here is large, and since the sum is converging, one can 
replace the summation by an integral which results in L17J 



Gel 



h 6E 



' Ec 



1 



+ 



1 



J\f-M* N* -N +1 



(25) 



Here we have included both the hole-like and electron-like contributions to the 
conductance and restored the explicit dependence on the gate voltage N; Eq. ( 1^51 
is valid for |7\A - TV* - l/2[ < 1/2 - 5E/Ec. Comparison of Eq. O with 
Eq. (EUl shows that the elastic co-tunneling dominates the electron transport al- 
ready at temperatures 

T<Tei = VEcSE, (26) 

which may exceed significantly the level spacing. 

Note that mesoscopic fluctuations of the elastic co-tunneling contribution to 
the conductance Gei are of the order of its average and get stronger when the gate 
voltage is tuned closer to the middle of the Coulomb blockade valley IT9l . Thus, 
although Gel is always positive, see Eq. d25l . the sample-specific value of Gei for 
a given gate voltage may vanish. 
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5. Kondo regime in transport through a blockaded quantum dot 

Among the Ec\M —N'*\/ 5E virtual states participating in the elastic co-tunneling 
through a blockaded dot, the top-most occupied single-particle level is special. If 
the number of electrons on the dot is odd, this level is filled by a single electron 
and is spin-degenerate. Therefore the ground state of the dot is characterized not 
only by the occupations of the single-particle energy levels, but also by the dot's 
spin. This opens a possibility of a co-tunneling process in which a transfer of 
an electron between the leads is accompanied by a flip of electron's spin with 
simultaneous flip of the spin on the dot, see Fig. Etc). The amplitude of such 
a process, calculated in the fourth order in tunneling matrix elements, diverges 
logarithmically when the energy uj of an incoming electron approaches 0. Since 
u; ~ T, the logarithmic singularity in the transmission amplitude translates into a 
dramatic enhancement of the conductance G across the dot at low temperatures: G 
may reach values as high as the quantum limit 2e^//i |20 21 1. This conductance 
enhancement is not really a surprise. Indeed, in the spin-flip co-tunneling process 
a quantum dot with odd N behaves as 5 = 1/2 magnetic impurity embedded into 
a tunneling barrier separating two massive conductors. It is known |22 1 since mid- 
60's that the presence of such impurities leads to zero-bias anomalies in tunneling 
conductance [23J, which are adequately explained 1241 in the context of the Kondo 
effect 0. 

For simplicity, we will assume here that the gate voltage J\f is tuned to the 
middle of the Coulomb blockade valley with N = odd electrons on the dot. The 
tunneling ^ mixes this state with states having A it 1 electrons. The electrostatic 
energies of these states are high (~ Ec), hence the transitions A — > A it 1 
are virtual, and can be taken into account perturbatively in the second order in 
tunneling amplitudes. The resulting effective Hamiltonian, valid at energies well 
below the threshold (~ 5E) for spin excitations within the dot, has the form 



and P is a projector onto the spin-degenerate ground state of an isolated dot. In 
the derivation of Eq. d27b we have replaced the tunneling amplitudes t^kn by their 
values tan at the Fermi level. In addition, we have dropped the potential scattering 
terms associated with usual elastic co-tunneling. The latter approximation is well 
justified when the conductances of the dot-leads junctions are small, Gq, <C /h, 
in which case G^i is also very small, see Eq. (I25t . 




Here 
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By SU (2) symmetry, the operators S„„' for any n and n' must be proportional 
to S = VSV with S being the operator of the total spin on the dot introduced in 
Eq. Q above: S„n' = A„„'S. Substitution into Eq. d27l then yields l25l 



with 



aks 



4 



Jaa' — tgn^nn'tgin' ■ ^^^^ 

nn' 

Under the conditions of applicability of the Constant Interaction Model Eq. (ISj, 
the effective exchange Hamiltonian d28l can be simplified even further. Indeed, in 
this case the ground state of an isolated dot is a singlet or a doublet, depending 
only on the parity of N. If N is odd, then S is spin- 1/2 operator. The entire spin 
of the dot is now due to the only singly occupied single-particle energy level in it 
(denoted by n = hereafter). The matrix A^n' then reduces to \nn' = ^nn'^no^ 
and Eq. (l29l yields 

4 

Jaa' = T^^aO^o'O • (30) 

The 2x2 Hermitian matrix J with elements given by Eq. ( BUI has an important 
property: since det J = 0, one of its eigenvalues vanishes, while the remaining 
eigenvalue, 

J = TrJ = -i-(|tio| + |t|ol), (31) 

is strictly positive. By an appropriate rotation in the R — L space the Hamiltonian 
(EHl can then be brought into the "block-diagonal" form 



H,ff = H[^] + H[^l;] (32) 



with 



ks 

Hm = Y.^k4s^ks + Ji^i'-s)- (34) 



ks 



Here = Ylkk'ss' ^Is (""ss' Z^) ''Pk's" operators ip and ip are certain Unear 

combinations of the original operators cj^^l- 

To get an idea about the physics of the Kondo model, let us first replace the 
operator in Eq. ( H^ by spin-1/2 operator S^. The ground state of the resulting 
Hamiltonian of two spins H = J(S^ • S) is obviously a singlet. The excited 
state (a triplet) is separated from the ground state by the energy gap J > 0. This 
separation can be interpreted as the binding energy of the singlet. Unlike in 



Kondo.tex; 2/02/2008; 0:31; p. 12 



COULOMB BLOCKADE AND KONDO EFFECT 



13 



this simple example, the operator in Eq. d34b is merely a spin density of the 
conduction electrons at the site of the "magnetic impurity". Because conduction 
electrons are freely moving in space, it is hard for the impurity to "capture" an 
electron and form a singlet. Yet, even a weak local exchange interaction suffices 
to form a singlet ground state L28,,22J- However, the binding energy of this singlet 
is given not by the exchange constant J, but by the so-called Kondo temperature 

Tk ^ 5Eexp{-l/i^J). (35) 

With the help of Eqs. dSTT i and (fTTl one obtains from Eq. d35l the estimate 

mf^U^^ ^. (36) 

'""{Tk) vJ Gl + Gr5E 

Since Ga ^ e'^/h, and Ec ^ SE, the r.h.s. of Eq. d36l is a product of two large 
parameters. Therefore, the Kondo temperature Tk is small compared to the mean 
level spacing in the dot: 

Tk < dE. (37) 

It is this separation of the energy scales that justifies the use of the effective low- 
energy Hamiltonian d27l for the description of the Kondo effect in a quantum dot 
system. The inequality ( 071 remains valid even if the conductances of the dot-leads 
junctions Ga are of the order of 2e^//i. However, in this case neither Eq. d3n i nor 
the estimate Eq. are applicable |30|. 

As it follows from Eqs. d33l and d34t . one of the "channels" (y?) of conduction 
electrons completely decouples^ from the dot, while the V^-particles are described 
by the standard single-channel antiferromagnetic Kondo model |5 1. Therefore, the 
thermodynamic properties of a quantum dot in the Kondo regime are identical to 
those of the conventional Kondo problem for a single magnetic impurity in a bulk 
metal; thermodynamics of the latter model is fully studied |26|. However, all the 
experiments addressing the Kondo effect in quantum dots test their transport prop- 
erties rather then thermodynamics. The electron current operator is not diagonal 
in the {xp, (/)) representation, and the contributions of these two sub-systems to the 
conductance are not additive. In the following we establish the relations of linear 
conductance and, in some special case, of the non-linear differential conductance 
as well, to the t-matrix of the conventional Kondo problem. 



^ It should be emphasized that the decoupling is characteristic of the Constant Interaction Model 
rather than a generic property. In general, both electronic channels are coupled to the dot. This, 
however, does not lead to qualitative changes in the results if spin of the dot is 1/2, see L25J . 
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5.1. LINEAR RESPONSE 

The linear conductance can be calculated from the Kubo formula 

p2 1 poo 

G=lim--/ dte-*([j(t),j(0)]), (38) 
where the particle current operator j is given by 

Here = X^fcs ^^aks^aks operator of the total number of electrons in the 
lead a. In order to take the full advantage of the decomposition (B^-d^. we 
need to rewrite j in terms of the operators ip, ip. These operators are related to the 
original operators cr^l representing the electrons in the right and left leads via 

\Vks J V ^Lks J 

where [/ is 2 x 2 unitary matrix that diagonalizes the matrix of the exchange 
constants d29l : U JW = diag{ J, 0}. The matrix U can be parametrized as 

(f = e'^oTy ^i<t>or, ^ (41) 

where are the PauU matrices acting in the R-L space (r+ = Tx + iTy transforms 
L to R) and the angle 6q satisfies tan^o = I^Lo/^RoI- It follows from Eqs. (l40l 
and (l4ni that, independently of (pQ, 

Nr-Nl = cos(20o) {n^ - N^) - sin(20o) {"t^l^ks + Asi^ks) ■ (42) 

ks 

Note that the current operator j entering the Kubo formula (l38l is to be calculated 
with the equilibrium Hamiltonian (l32t - (l34t . Since both and iV^ commute 
with Hf,ff, the first term in Eq. d42t makes no contribution to j. When the second 
term in d?2t is substituted into Eq. (l39t and then into Eq. d38l . the result, after 
integration by parts, can be expressed via 2-particle correlation functions such as 
(^ip^ {t)ip{t)(p^ {0)'ijj{0)y Due to the block-diagonal structure of Hf,ff, see Eq. (l32t . 
these correlation functions factorize into products of single-particle correlation 
functions describing the (free) (^-particles and the (interacting) V'-particles. The 
result of the evaluation of the Kubo formula can then be written as^ 

G = Go f dio i-df/dio) ^ [-^^I™ ^-(^)] • (43) 

s 

^ Further details about this calculation can be found, e.g., in the Appendix B of Ref. 1271 . 
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Here 



— W{2e,) = — 



2 



Go 




(44) 



m 



f{uj) is the Fermi function, and Tg is the t-matrix for the Kondo model Eq. d34t . 
The t-matrix is related to the exact retarded Green function of the i/'-particles in 
the usual way: 



Here G1{uj) = {lo — + iO) ^ and Gks,k's'{^) is the Fourier transform of 



where (...) stands for the thermodynamic averaging with the Hamiltonian d34l i. 
In writing Eq. (l45l we took into account the conservation of the total spin (which 
implies that Gks,k's' is diagonal in s, s') and the fact that the interaction in Eq. d34l i 
is local (which in turn means that the t-matrix is independent of k and k'). 

5.1.1. Weak coupling regime: Tk -^T <^ 5E 

Kondo effect becomes important at low temperatures T Tk ^ 6E. The ex- 
change term in the Hamiltonian d34b . however, describes transitions between the 
electronic states within the band of the width 2Dq with Dq ~ 6E centered at the 
Fermi level ep- The transitions between the states close to ei? and the states within 
a narrow strip of energies of the width 6D near the edges of the band are associated 
with high energy deficit ~ 6E 3> T. Hence, these transitions are virtual and their 
influence on the states near ep can be taken into account perturbatively in the 
second order. This yields an effective Hamiltonian acting within the band of a 
reduced width D = Dq — 6D, which turns out to have the same form as Eq. d34l i. 
but with a modified value of the exchange amplitude J |32|. Successive reduc- 
tions of the bandwidth by small steps 6D can be viewed as a continuous process 
during which the initial Hamiltonian Eq. d34T) is transformed to an effective low- 
energy Hamiltonian acting within the band of the width D ^ Dq. The evolution 
of the exchange amplitude during this transformation (known as the poor man's 
scaling |32l) can be cast into the form of an equation 



The renormalization described by this equation can be continued until the band- 
width D becomes of the order of the temperature T or relevant electron energy 
|a;|. The resulting effective exchange constant is lo- and T-dependent Ii32ii33i . 



Gks,k'Aio) = Gl{u:) + Gl{u:) Gl 



(45) 



Gks,k's'it) = -ie{t) {{^ksit)A's]) 



— = vj\ C = \u{PqID). 



max{|w[, r} 



n -1 



vJ{uj, T) = In 



Tk 



(46) 
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Calculation of the t-matrix in the lowest (second) order in uJ{uj, T) results in 

- 7ri.Imr,(a;) = — [vJ{uJ,T)f . (47) 

Substitution of (l47l and d46l into Eq. (l43l and evaluation of the integral over ui 
with the logarithmic accuracy then yield for the conductance across the dot 

Corrections to Eq. (l48l contain higher powers of 1/ h\{T /Tk)- Thus, it is often 
said that Eq. d48l represents the conductance in the leading logarithmic approxi- 
mation. 



G = Go ^^ 1 2 ' Tk-^T-^SE. (48) 



5.1.2. Strong coupling regime: T <^ Tk 

As temperature approaches T^-, the leading logarithmic approximation result di- 
verges, see Eq. (l48l . This divergence signals the failure of the approximation, as 
the conductance in any case can not grow higher than 2e^ /h. To obtain a more 
precise bound, we consider in this section the conductance in the strong coupling 
regime T «; T^. 

We start with the zero-temperature limit T = 0. Since, as discussed above, 
the ground state of the Hamiltonian d34l i is not degenerate, the scattering of con- 
duction electrons by a quantum dot is completely characterized by the scattering 
phase shifts 5s for electrons with spin s at the Fermi level. The t-matrix is then 
given by the standard scattering theory expression 

- vrz^TJO) = - ie^'^^ - l) = e'^^ sin 5s. (49) 
2i \ J 

In order to find the two phase shifts 5s, we need two independent relations. 
The first one follows from the invariance of the Kondo Hamiltonian (1^ under 
the particle-hole transformation — > sV'lfc _s (here s = ±1 for spin-up/down 
electrons and k is counted from kp). The particle-hole symmetry implies the 
relation for the t-matrix 

r,(u;) = -r:,(-^), 

valid at all lo and T. In view of Eq. (l49l . it translates into the relation for the phase 
shifts at the Fermi level (a; = 0): 

5^ + = (modvr). (50) 

Note that, as it is obvious from Eq. ( l49l . the phase shifts are defined only modulo 
TT (that is, 5s and (5s + vr describe equivalent scattering states). This ambiguity can 
be removed by setting the values of the phase shifts corresponding to J = in 
Eq. (l34l i to zero. With this convention, Eq. (l50l becomes 

J^ + 5^ = 0. (51) 
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The second relation follows from the fact that the ground state of our system 
is a singlet at J > 0. In the absence of exchange (J = 0) and at T = 0, an 
infinitesimally small magnetic field acting on the dot's spin, 

5H[^Ij] = BS'\ B +0, 

would polarize it. Since the free electron gas has zero spin in the ground state, the 
total spin in any large but finite region of space V surrounding the dot coincides 
with the spin of the dot, {S^)j=o = —1/2. If the exchange with electron gas is 
now turned on, J > 0, the infinitesimally weak field will not prevent the formation 
of a singlet ground state. In this state, the total spin within the region V is zero. 
Such change of the spin is possible only if the numbers of spin-up and spin-down 
electrons in this region have changed in order to compensate for the spin of the 
dot: 

6N^ - 5Ni = 1, 6Ns = (iV,)j>o - (/^.)j=o. (52) 

Here Ns are operators of the number of electrons with spin s within the region 
V. By the Friedel sum rule, 6Ns are related to the scattering phase shifts: 6Ns = 
6s /tt. Eq. d52b then gives 

6^-Si= IT. (53) 

Combining Eqs. dSTT i and d53l . we find \5s\ = 7i"/2. Eqs. d43t and d49b then yield 
for zero-temperature conductance across the dot |20| 

G'(0) = Go^J]sin2 5, = Go. (54) 

s 

Thus, the grows of the conductance across the dot with lowering the temperature 
is limited only by the value of Gq. This value, see Eq. (l44T i. depends only on the 
ratio of the tunneling amplitudes |tLo/iiJo| between the leads and the last occupied 
single-particle energy level in the dot. If |tLo| = I^RoL the conductance at T = 
will reach the maximal value G = 2e^//i allowed by quantum mechanics |20|. 

The ground state of the Kondo Hamiltonian (l34l is a singlet formed by the 
spin of the dot and a spin made up of the spins of conduction electrons. Finite- 
temperature correction to Eq. d55l can be found by considering virtual transitions 
from the singlet ground state to excited states in which the singlet is broken 
up l34l . The transitions can be studied by an expansion in inverse powers of the 
singlet binding energy Tk- The reader is refen^ed to the original papers ['341 for 
the details about this approach. Here we just quote the result for the imaginary 
part of the t-matrix ll35l : 



^z.Imr,(c^) = 1 —2 , \uj\,T^Tk (55) 



Substitution of Eq. d55l into Eq. (l43l yields 



G = Gn 



{ttT/Tk) 



T^Tk (56) 
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Accordingly, corrections to the conductance are quadratic in temperature - a typ- 
ical result for the Fermi liquid theory |34|. 

The weak-coupling (T ^ Tk) and the strong-coupling (T ^ Tk) asymptotes 
of the conductance have a strikingly different structure. Nevertheless, since the 

Kondo effect is a crossover phenomenon rather than a phase transition L26 28..29J . 

the dependence G{T) is a smooth and featureless |36| function throughout the 
crossover region T Tk- 

Finally, note that both Eqs. (l48l and d56t have been obtained here for the 
particle-hole symmetric model d34l . This approximation is equivalent to neglect- 
ing the elastic co-tunneling contribution to the conductance Gei- The asymptotes 
(l48l . (l54l . however, remain qualitatively correct 1*251 as long as Gei/Go <^ 1. 
The overall temperature dependence of the linear conductance in the Coulomb 
blockade valley in the presence of the Kondo effect is sketched in Fig. 3. 



G 

Gq 



(56)' 



Ga 



Ge 





(48) 


(25) 


(20) 


(17) 


(12) 

























T 



K 



5E 



Tel 



Ec 



T 



Figure 3. Schematic temperature dependence of the conductance in the middle of the Coulomb 
blockade valley with A'^ — odd electrons on the dot. The numbers in brackets refer to the 
corresponding equations in the text. 



5.2. BEYOND LINEAR RESPONSE 

In order to study transport through a quantum dot away from equilibrium we add 
to the effective Hamiltonian d32b - d34b a term 

Hv = ^ [Nl - Nr) 

describing a finite voltage bias V applied between the left and right electrodes. 
Here we will evaluate the current across the dot at arbitrary V but under the simpli- 
fying assumption that the dot-lead junctions are strongly asymmetric: Gl <SC Gr. 
Under this condition the angle in Eq. d^Tl is small, 

Oo ~ \tLo/tm\ <■ 1- 
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Expanding Eq. (l42l to the linear order in 6*0 we obtain 

Hv{Oo) = ^{n^- N^) + eVeoA, (57) 

ks 

The first term in the r.h.s. of Eq. dSTl describes the voltage bias between the 
reservoirs of ip and ^/^-particles, while the second term has an appearance of the 
^-conserving "tunneling" with ^-dependent "tunneling amplitude". These two 
contributions correspond to the two terms in the r.h.s. of Eq. d42b . Similarly, the 
current operator j, see Eq. ( l39t . can be also split into two contributions: 

(58) 

= ieyeo^'^lA, + H.c., 

ks 

It turns out that 6j in Eq. ( I5HI makes no contribution to the average current in the 
leading (second) order in ^o- Indeed, in this order 

where (...)o denotes thermodynamic averaging with the Hamiltonian Hq = Hf,ff+ 
Hv{0)- The thermodynamic (equilibrium) averaging is well defined here because 
the "tunneling" term in Eq. (1571 is absent at 6*0 = 0. Taking into account that 

([A{t'),A{t)fj^ depends only ont - t' = t, we find 

{Sj) = ^OleV- dr ([A(r), ^(0)])^ ^ 0. 

To the contrary, the term jo in Eq. d58l makes a finite contribution to the average 
current. Evaluation of this contribution in the leading order in proceeds simi- 
larly to the calculation of the tunneling current between two massive conductors 
(see, e.g., IIbTI '). and yields for the differential conductance across the dot 

dl/dV = Go 1 duo {-df/duj) ]- ^ [-vrz/Im r,(w + eV)] (59) 

J ^ s 

with 



J 

jo 
53 



h + oj, 
i 



2 L" 



'4- 
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Note that the V Hmit of Eq. d59t coincides with the small 6*0 limit of the 
linear response result Eq. Using now Eqs. (l?^.(H7l.(l55l. and ^% we obtain 
for the differential conductance atT ^ eV <^ 5E 



dl_ 

dV 



Go 



1 



3 /eV 
2 IfV 



37rVl6 
MeV/TK)f 



eV < Tk 



eV->TK 



(60) 



Thus, a large voltage bias eV ^ T has qualitatively the same destructive effect 
on the Kondo physics as the temperature does. If temperature T exceeds the 
bias, T » eV, the differential conductance dl/dV coincides with the linear 
conductance G, see Eqs. (I48l . d56t above. 



6. Conclusion 



In the simplest form of the Kondo effect considered in this article, a quantum dot 
behaves essentially as a magnetic impurity with spin 1/2 coupled via exchange 
interaction to two conducting leads |24|. However, the characteristic energy scale 
for the intra-dot excitations is much smaller than the corresponding scale for a real 
magnetic impurity. This allows one to induce some degeneracies in the ground 
state of a dot which are more exotic than just the spin degeneracy considered 
above. One of the possibilities is to create a degeneracy between a singlet state 
and a component of a triplet state by applying a magnetic field to a dot with an 
even number of electrons. Depending on the relation between the Zeeman energy 
and the value of the diamagnetic shift of orbital levels in a magnetic field, this 
degeneracy gives rise to various "flavors" of the Kondo effect. A review of the 
corresponding experiments and theory can be found in 1 37 1 

In our discussion of the Kondo effect we assumed that the dot remains close to 
an equilibrium state even under the conditions of an electron transport experiment. 
In fact, this limitation led us to consider a very asymmetric setup of a quantum 
dot device in Sec. 15. 21 where the non-linear electron transport is discussed. One of 
the advantages of quantum dots, however, is that these devices allow one to study 
the Kondo effect in truly out-of-equilibrium conditions. Such conditions may be 
created by applying a significant dc bias between the leads, or by irradiating a 
dot by microwaves |38|. In the latter case one can monitor the "health" of the 
Kondo effect by measuring the linear dc conductance at the same time. It turns 
out that microwaves suppress the effect |38| by destroying the coherence of the 
spin state of the dot |39|. In the former case, one can apply a strong magnetic 
field in addition to biasing the dot. Zeeman splitting then leads to peaks in the 
differential conductance at a finite bias [40"4T1. Interest to this problem was 
recently revived, see, e.g., [42 1 and references therein. 
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These are only a few out of many possible extensions of the simple model 
discussed in this review. 
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